Virulent Newcastle disease virus genotypes V.3, VII.2, and XIII.1.1 and their coinfections with infectious bronchitis viruses and other avian pathogens in backyard chickens in Tanzania.

Oropharyngeal (OP) and cloacal (CL) swabs from 2049 adult backyard chickens collected at 12 live bird markets, two each in Arusha, Dar es Salaam, Iringa, Mbeya, Morogoro and Tanga regions of Tanzania were screened for Newcastle disease virus (NDV) using reverse transcription real-time PCR (rRT-PCR). The virus was confirmed in 25.23% of the birds (n = 517; rRT-PCR CT ≤ 30), with the highest positivity rates observed in birds from Dar es Salaam region with higher prevalence during the dry season (September-November 2018) compared to the rainy season (January and April-May 2019). Next-generation sequencing of OP/CL samples of 20 out of 32 birds that had high amounts of viral RNAs (CT ≤ 25) resulted in the assembly of 18 complete and two partial genome sequences (15,192 bp and 15,045-15,190 bp in length, respectively) of NDV sub-genotypes V.3, VII.2 and XIII.1.1 (n = 1, 13 and 4 strains, respectively). Two birds had mixed NDV infections (V.3/VII.2 and VII.2/XIII.1.1), and nine were coinfected with viruses of families Astroviridae, Coronaviridae, Orthomyxoviridae, Picornaviridae, Pneumoviridae, and Reoviridae. Of the coinfecting viruses, complete genome sequences of two avastroviruses (a recombinant chicken astrovirus antigenic group-Aii and avian nephritis virus genogroup-5) and two infectious bronchitis viruses (a turkey coronavirus-like recombinant and a GI-19 virus) were determined. The fusion (F) protein F1/F2 cleavage sites of the Tanzanian NDVs have the consensus motifs 112 RRRKR↓F 117 (VII.2 strains) and 112 RRQKR↓F 117 (V.3 and XIII.1.1 strains) consistent with virulent virus; virulence was confirmed by intracerebral pathogenicity index scores of 1.66-1.88 in 1-day-old chicks using nine of the 20 isolates. Phylogenetically, the complete F-gene and full genome sequences regionally cluster the Tanzanian NDVs with, but distinctly from, other strains previously reported in eastern and southern African countries. These data contribute to the understanding of NDV epidemiology in Tanzania and the region.

Molecular Characterization of Complete Genome Sequence of an Avian Coronavirus Identified in a Backyard Chicken from Tanzania.

A complete genome sequence of an avian coronavirus (AvCoV; 27,663 bp excluding 3' poly(A) tail) was determined using nontargeted next-generation sequencing (NGS) of an oropharyngeal swab from a backyard chicken in a live bird market in Arusha, Tanzania. The open reading frames (ORFs) of the Tanzanian strain TZ/CA127/19 are organized as typical of gammaCoVs (Coronaviridae family): 5'UTR-[ORFs 1a/1b encoding replicase complex (Rep1ab) non-structural peptides nsp2-16]-[spike (S) protein]-[ORFs 3a/3b]-[small envelop (E) protein]-[membrane (M) protein]-[ORFs 4a/4c]-[ORFs 5a/5b]-[nucleocapsid (N) protein]-[ORF6b]-3'UTR. The structural (S, E, M and N) and Rep1ab proteins of TZ/CA127/19 contain features typically conserved in AvCoVs, including the cleavage sites and functional motifs in Rep1ab and S. Its genome backbone (non-spike region) is closest to Asian GI-7 and GI-19 infectious bronchitis viruses (IBVs) with 87.2-89.7% nucleotide (nt) identities, but it has a S gene closest (98.9% nt identity) to the recombinant strain ck/CN/ahysx-1/16. Its 3a, 3b E and 4c sequences are closest to the duck CoV strain DK/GD/27/14 at 99.43%, 100%, 99.65% and 99.38% nt identities, respectively. Whereas its S gene phylogenetically cluster with North American TCoVs and French guineafowl COVs, all other viral genes group monophyletically with Eurasian GI-7/GI-19 IBVs and Chinese recombinant AvCoVs. Detection of a 4445 nt-long recombinant fragment with breakpoints at positions 19,961 and 24,405 (C- and N-terminus of nsp16 and E, respectively) strongly suggested that TZ/CA127/19 acquired its genome backbone from an LX4-type (GI-19) field strain via recombination with an unknown AvCoV. This is the first report of AvCoV in Tanzania and leaves unanswered the questions of its emergence and the biological significance.

Genome Sequences and Characterization of Chicken Astrovirus and Avian Nephritis Virus from Tanzanian Live Bird Markets.

The enteric chicken astrovirus (CAstV) and avian nephritis virus (ANV) are the type species of the genus Avastrovirus (AAstV; Astroviridae family), capable of causing considerable production losses in poultry. Using next-generation sequencing of a cloacal swab from a backyard chicken in Tanzania, we assembled genome sequences of ANV and CAstV (6918 nt and 7318 nt in length, respectively, excluding poly(A) tails, which have a typical AAstV genome architecture (5'-UTR-ORF1a-ORF1b-ORF2-'3-UTR). They are most similar to strains ck/ANV/BR/RS/6R/15 (82.72%) and ck/CAstV/PL/G059/14 (82.23%), respectively. Phylogenetic and sequence analyses of the genomes and the three open reading frames (ORFs) grouped the Tanzanian ANV and CAstV strains with Eurasian ANV-5 and CAstV-Aii viruses, respectively. Compared to other AAstVs, the Tanzanian strains have numerous amino acid variations (substitutions, insertions and deletions) in the spike region of the capsid protein. Furthermore, CAstV-A has a 4018 nt recombinant fragment in the ORF1a/1b genomic region, predicted to be from Eurasian CAstV-Bi and Bvi parental strains. These data should inform future epidemiological studies and options for AAstV diagnostics and vaccines.

Viral proteogenomic and expression profiling during productive replication of a skin-tropic herpesvirus in the natural host.

Efficient transmission of herpesviruses is essential for dissemination in host populations; however, little is known about the viral genes that mediate transmission, mostly due to a lack of natural virus-host model systems. Marek's disease is a devastating herpesviral disease of chickens caused by Marek's disease virus (MDV) and an excellent natural model to study skin-tropic herpesviruses and transmission. Like varicella zoster virus that causes chicken pox in humans, the only site where infectious cell-free MD virions are efficiently produced is in epithelial skin cells, a requirement for host-to-host transmission. Here, we enriched for heavily infected feather follicle epithelial skin cells of live chickens to measure both viral transcription and protein expression using combined short- and long-read RNA sequencing and LC/MS-MS bottom-up proteomics. Enrichment produced a previously unseen breadth and depth of viral peptide sequencing. We confirmed protein translation for 84 viral genes at high confidence (1% FDR) and correlated relative protein abundance with RNA expression levels. Using a proteogenomic approach, we confirmed translation of most well-characterized spliced viral transcripts and identified a novel, abundant isoform of the 14 kDa transcript family via IsoSeq transcripts, short-read intron-spanning sequencing reads, and a high-quality junction-spanning peptide identification. We identified peptides representing alternative start codon usage in several genes and putative novel microORFs at the 5' ends of two core herpesviral genes, pUL47 and ICP4, along with strong evidence of independent transcription and translation of the capsid scaffold protein pUL26.5. Using a natural animal host model system to examine viral gene expression provides a robust, efficient, and meaningful way of validating results gathered from cell culture systems.

Complete Genome Sequences of Avian Metapneumovirus Subtype B Vaccine Strains from Brazil.

Avian metapneumovirus (aMPV) causes a highly contagious upper respiratory and reproductive disease in chickens, turkeys, and ducks. Here, complete genome sequences of aMPV-B vaccine strains BR/1890/E1/19 (PL21, Nemovac; Boehringer Ingelheim Animal Health, Brazil) and BR/1891/E2/19 (1062; Hipraviar, France) were sequenced and compared with the pathogenic field strain VCO3/60616.

Unique Variants of Avian Coronaviruses from Indigenous Chickens in Kenya.

The avian gamma-coronavirus infectious bronchitis virus (AvCoV, IBV; Coronaviridae family) causes upper respiratory disease associated with severe economic losses in the poultry industry worldwide. Here, we report for the first time in Kenya and the Eastern African region two novel AvCoVs, designated IBV/ck/KE/1920/A374/2017 (A374/17) and AvCoV/ck/KE/1922/A376/2017 (A376/17), inadvertently discovered using random nontargeted next-generation sequencing (NGS) of cloacal swabs collected from indigenous chickens. Despite having genome organization (5'UTR-[Rep1a/1ab-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR), canonical conservation of essential genes and size (~27.6 kb) typical of IBVs, the Kenyan isolates do not phylogenetically cluster with any genotypes of the 37 IBV lineages and 26 unique variants (UVs). Excluding the spike gene, genome sequences of A374/17 and A376/17 are only 93.1% similar to each other and 86.7-91.4% identical to genomes of other AvCoVs. All five non-spike genes of the two isolates phylogenetically cluster together and distinctly from other IBVs and turkey coronaviruses (TCoVs), including the indigenous African GI-26 viruses, suggesting a common origin of the genome backbone of the Kenyan isolates. However, isolate A376/17 contains a TCoV-like spike (S) protein coding sequence and is most similar to Asian TCoVs (84.5-85.1%) compared to other TCoVs (75.6-78.5%), whereas isolate A374/17 contains an S1 gene sequence most similar to the globally distributed lineage GI-16 (78.4-79.5%) and the Middle Eastern lineage GI-23 (79.8-80.2%) viruses. Unanswered questions include the actual origin of the Kenyan AvCoVs, the potential pathobiological significance of their genetic variations, whether they have indeed established themselves as independent variants and subsequently spread within Kenya and to the neighboring east/central African countries that have porous live poultry trade borders, and whether the live-attenuated Mass-type (lineage GI-1)-based vaccines currently used in Kenya and most of the African countries provide protection against these genetically divergent field variants.

Comparable outcomes from long and short read random sequencing of total RNA for detection of pathogens in chicken respiratory samples.

Co-infections of avian species with different RNA viruses and pathogenic bacteria are often misdiagnosed or incompletely characterized using targeted diagnostic methods, which could affect the accurate management of clinical disease. A non-targeted sequencing approach with rapid and precise characterization of pathogens should help respiratory disease management by providing a comprehensive view of the causes of disease. Long-read portable sequencers have significant potential advantages over established short-read sequencers due to portability, speed, and lower cost. The applicability of short reads random sequencing for direct detection of pathogens in clinical poultry samples has been previously demonstrated. Here we demonstrate the feasibility of long read random sequencing approaches to identify disease agents in clinical samples. Experimental oropharyngeal swab samples (n = 12) from chickens infected with infectious bronchitis virus (IBV), avian influenza virus (AIV) and Mycoplasma synoviae (MS) and field-collected clinical oropharyngeal swab samples (n = 11) from Kenyan live bird markets previously testing positive for Newcastle disease virus (NDV) were randomly sequenced on the MinION platform and results validated by comparing to real time PCR and short read random sequencing in the Illumina MiSeq platform. In the swabs from experimental infections, each of three agents in every RT-qPCR-positive sample (Ct range 19-34) was detectable within 1 h on the MinION platform, except for AIV one agent in one sample (Ct = 36.21). Nine of 12 IBV-positive samples were assigned genotypes within 1 h, as were five of 11 AIV-positive samples. MinION relative abundances of the test agent (AIV, IBV and MS) were highly correlated with RT-qPCR Ct values (R range-0.82 to-0.98). In field-collected clinical swab samples, NDV (Ct range 12-37) was detected in all eleven samples within 1 h of MinION sequencing, with 10 of 11 samples accurately genotyped within 1 h. All NDV-positive field samples were found to be co-infected with one or more additional respiratory agents. These results demonstrate that MinION sequencing can provide rapid, and sensitive non-targeted detection and genetic characterization of co-existing respiratory pathogens in clinical samples with similar performance to the Illumina MiSeq.

Detection and Genome Sequence Analysis of Avian Metapneumovirus Subtype A Viruses Circulating in Commercial Chicken Flocks in Mexico.

Avian metapneumoviruses (aMPV subtypes A-D) are respiratory and reproductive pathogens of poultry. Since aMPV-A was initially reported in Mexico in 2014, there have been no additional reports of its detection in the country. Using nontargeted next-generation sequencing (NGS) of FTA card-spotted respiratory samples from commercial chickens in Mexico, seven full genome sequences of aMPV-A (lengths of 13,288-13,381 nucleotides) were de novo assembled. Additionally, complete coding sequences of genes N (n = 2), P and M (n = 7 each), F and L (n = 1 each), M2 (n = 6), SH (n = 5) and G (n = 2) were reference-based assembled from another seven samples. The Mexican isolates phylogenetically group with, but in a distinct clade separate from, other aMPV-A strains. The genome and G-gene nt sequences of the Mexican aMPVs are closest to strain UK/8544/06 (97.22-97.47% and 95.07-95.83%, respectively). Various amino acid variations distinguish the Mexican isolates from each other, and other aMPV-A strains, most of which are in the G (n = 38), F (n = 12), and L (n = 19) proteins. Using our sequence data and publicly available aMPV-A data, we revised a previously published rRT-PCR test, which resulted in different cycling and amplification conditions for aMPV-A to make it more compatible with other commonly used rRT-PCR diagnostic cycling conditions. This is the first comprehensive sequence analysis of aMPVs in Mexico and demonstrates the value of nontargeted NGS to identify pathogens where targeted virus surveillance is likely not routinely performed.

Genome Sequence Variations of Infectious Bronchitis Virus Serotypes From Commercial Chickens in Mexico.

New variants of infectious bronchitis viruses (IBVs; Coronaviridae) continuously emerge despite routine vaccinations. Here, we report genome sequence variations of IBVs identified by random non-targeted next generation sequencing (NGS) of vaccine and field samples collected on FTA cards from commercial flocks in Mexico in 2019-2021. Paired-ended sequencing libraries prepared from rRNA-depleted RNAs were sequenced using Illumina MiSeq. IBV RNA was detected in 60.07% (n = 167) of the analyzed samples, from which 33 complete genome sequences were de novo assembled. The genomes are organized as 5'UTR-[Rep1a-Rep1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR, except in eight sequences lacking non-structural protein genes (accessory genes) 4b, 4c, and 6b. Seventeen sequences have auxiliary S2' cleavage site located 153 residues downstream the canonically conserved primary furin-specific S1/S2 cleavage site. The sequences distinctly cluster into lineages GI-1 (Mass-type; n = 8), GI-3 (Holte/Iowa-97; n = 2), GI-9 (Arkansas-like; n = 8), GI-13 (793B; n = 14), and GI-17 (California variant; CAV; n = 1), with regional distribution in Mexico; this is the first report of the presence of 793B- and CAV-like strains in the country. Various point mutations, substitutions, insertions and deletions are present in the S1 hypervariable regions (HVRs I-III) across all 5 lineages, including in residues 38, 43, 56, 63, 66, and 69 that are critical in viral attachment to respiratory tract tissues. Nine intra-/inter-lineage recombination events are present in the S proteins of three Mass-type sequences, two each of Holte/Iowa-97 and Ark-like sequence, and one each of 793B-like and CAV-like sequences. This study demonstrates the feasibility of FTA cards as an attractive, adoptable low-cost sampling option for untargeted discovery of avian viral agents in field-collected clinical samples. Collectively, our data points to co-circulation of multiple distinct IBVs in Mexican commercial flocks, underscoring the need for active surveillance and a review of IBV vaccines currently used in Mexico and the larger Latin America region.

Surveillance and Genetic Characterization of Virulent Newcastle Disease Virus Subgenotype V.3 in Indigenous Chickens from Backyard Poultry Farms and Live Bird Markets in Kenya.

Kenyan poultry consists of ~80% free-range indigenous chickens kept in small flocks (~30 birds) on backyard poultry farms (BPFs) and they are traded via live bird markets (LBMs). Newcastle disease virus (NDV) was detected in samples collected from chickens, wild farm birds, and other domestic poultry species during a 2017-2018 survey conducted at 66 BPFs and 21 LBMs in nine Kenyan counties. NDV nucleic acids were detected by rRT-PCR L-test in 39.5% (641/1621) of 1621 analyzed samples, of which 9.67% (62/641) were NDV-positive by both the L-test and a fusion-test designed to identify the virulent virus, with a majority being at LBMs (64.5%; 40/62) compared to BPFs (25.5%; 22/62). Virus isolation and next-generation sequencing (NGS) on a subset of samples resulted in 32 complete NDV genome sequences with 95.8-100% nucleotide identities amongst themselves and 95.7-98.2% identity with other east African isolates from 2010-2016. These isolates were classified as a new sub-genotype, V.3, and shared 86.5-88.9% and 88.5-91.8% nucleotide identities with subgenotypes V.1 and V.2 viruses, respectively. The putative fusion protein cleavage site (113R-Q-K-R↓F 117) in all 32 isolates, and a 1.86 ICPI score of an isolate from a BPF chicken that had clinical signs consistent with Newcastle disease, confirmed the high virulence of the NDVs. Compared to genotypes V and VI viruses, the attachment (HN) protein of 18 of the 32 vNDVs had amino acid substitutions in the antigenic sites. A time-scaled phylogeographic analysis suggests a west-to-east dispersal of the NDVs via the live chicken trade, but the virus origins remain unconfirmed due to scarcity of continuous and systematic surveillance data. This study reveals the widespread prevalence of vNDVs in Kenyan backyard poultry, the central role of LBMs in the dispersal and possibly generation of new virus variants, and the need for robust molecular epidemiological surveillance in poultry and non-poultry avian species.

Runting and Stunting Syndrome in Broiler Chickens: Histopathology and Association With a Novel Picornavirus.

Runting stunting syndrome (RSS) in commercial chickens has been reported worldwide, and although several studies have attempted to clarify the cause and describe the lesions, there are gaps in knowledge of the epidemiology, pathogenesis, and etiology. The study objective was to use commercial chicks naturally affected by RSS to describe the histologic changes of RSS in all segments of the small intestine in chicks of different ages and to identify viral gene sequences in affected chicks and their association with histologic lesions. Chicks lacking clinical signs but from the same houses and from unaffected houses were used as controls. The average weight of affected chicks was significantly lower than expected for their flocks. Macroscopically, the small intestines had paler serosa, with watery, mucoid, or foamy contents and poorly digested food. Histologic lesions were characterized by necrotic crypts, crypt dilation, and flattening of the crypt epithelium. Histomorphometry of the intestines revealed villous atrophy especially in the jejunum and ileum. Histologic changes in other organs were not observed. Random next-generation sequencing of total RNA extracted from formalin-fixed paraffin-embedded tissues detected avian nephritis virus, avian rotavirus, and picornavirus in jejunal segments from 7-day-old chicks. No viruses were detected in the jejunum of 1-day-old chicks. Detection of picornaviral reads was significantly associated (P < .05) with histologic lesions of RSS. Sequence analysis of the picornavirus revealed genetic similarity with the genus Gallivirus. Using in situ hybridization for galliviral nucleic acid sequences, the signal was associated with crypt lesion severity, although signal was detected both in chicks with and without RSS.

Comparative Molecular Characterization of Three Gallid alphaherpesvirus Type 3 Strains 301B/1, HPRS24, and SB-1.

Marek's disease (MD) is a highly contagious lymphoproliferative disease of chickens caused by Gallid alphaherpesvirus type 2. Gallid alphaherpesvirus type 3 (GaHV-3) strain 301B/1 was previously shown to be an effective MD vaccine with synergistic efficacy when used as a bivalent vaccine with turkey herpesvirus. Since the nucleotide sequences of only two GaHV-3 strains have been determined, we sought to sequence the 301B/1 genome using Illumina MiSeq technology. Phylogenomic analysis indicated that 301B/1 is more closely related to other GaHV-3 strains (SB-1 and HPRS24) than to virulent or attenuated strains of GaHV-2. One hundred and twenty-six open reading frames (ORFs) have been identified within the 301B/1 genome with 108 ORFs showing a high degree of similarity to homologs found in the genomes of SB-1 and HPRS24; 14 ORFs are highly homologous (> 90% identity) with the corresponding ORFs within the SB-1 genome. The R-LORF8 and R-LORF9 genes are the most dissimilar to the collinear genes found in the SB-1 genome but are highly homologous (99%-100% identity) with those within the HPRS24 genome. Overall the 301B/1 genome is most similar to the SB-1 virus genome (99.1%) and to a lesser degree with the HPRS24 virus genome (97.7%). However, six 301B/1 ORFs (UL47, UL48, UL52, pp38, ICP4, and US10) have been identified that contain nonsynonymous substitutions relative to homologs found in the SB-1 genome. Notably, unlike the avian retrovirus long terminal repeat sequences found within the SB-1 genome, none were identified within the 301B/1 genome.

Caracterización molecular comparativa de cepas de Alfaherpesvirus del pollo tipo 3 cepas 301B/1, HPRS24 y SB-1. La enfermedad de Marek (MD) es una enfermedad linfoproliferativa altamente contagiosa de los pollos causada por el Alfaherpesvirus del pollo tipo 2. Se demostró previamente que la cepa 301B/1 del Alfaherpesvirus del pollo tipo 3 (GaHV-3) es una vacuna eficaz contra la enfermedad de Marek con eficacia sinérgica cuando se usa como una vacuna bivalente con el herpesvirus del pavo. Dado que se han determinado las secuencias de nucleótidos de solo dos cepas de GaHV-3, se buscó secuenciar el genoma de la cepa 301B/1 utilizando la tecnología Illumina MiSeq. El análisis filogenómico indicó que la cepa 301B/1 está más estrechamente relacionado con otras cepas de GaHV-3 (SB-1 y HPRS24) en comparación con cepas virulentas o atenuadas de GaHV-2. Se han identificado 126 marcos de lectura continuos (ORF) dentro del genoma de la cepa 301B/1 con 108 marcos de lectura continuos que muestran un alto grado de similitud con los secuencias homólogas encontrados en los genomas de las cepas SB-1 y HPRS24; 14 marcos de lectura continuo son altamente similares (> 90% de identidad) con los correspondientes dentro del genoma de SB-1. Los genes R-LORF8 y R-LORF9 fueron los más diferentes a los genes colineales encontrados en el genoma de SB-1, pero son altamente similares (99% -100% de identidad) con aquellos dentro del genoma HPRS24. En general, el genoma de la cepa 301B/1 es más similar al genoma del virus SB-1 (99.1%) y en menor grado con el genoma del virus HPRS24 (97.7%). Sin embargo, se han identificado seis marcos de lectura continuos en 301B/1 (UL47, UL48, UL52, pp38, ICP4 y US10) que contienen sustituciones no sinónimas en relación con las secuencias homólogas encontradas en el genoma SB-1. Notablemente, a diferencia de las secuencias repetidas terminales largas del retrovirus aviar encontradas dentro del genoma de SB-1, ninguna se identificó dentro del genoma 301B/1.

Rapid, multiplexed, whole genome and plasmid sequencing of foodborne pathogens using long-read nanopore technology.

U.S. public health agencies have employed next-generation sequencing (NGS) as a tool to quickly identify foodborne pathogens during outbreaks. Although established short-read NGS technologies are known to provide highly accurate data, long-read sequencing is still needed to resolve highly-repetitive genomic regions and genomic arrangement, and to close the sequences of bacterial chromosomes and plasmids. Here, we report the use of long-read nanopore sequencing to simultaneously sequence the entire chromosome and plasmid of Salmonella enterica subsp. enterica serovar Bareilly and Escherichia coli O157:H7. We developed a rapid and random sequencing approach coupled with de novo genome assembly within a customized data analysis workflow that uses publicly-available tools. In sequencing runs as short as four hours, using the MinION instrument, we obtained full-length genomes with an average identity of 99.87% for Salmonella Bareilly and 99.89% for E. coli in comparison to the respective MiSeq references. These nanopore-only assemblies provided readily available information on serotype, virulence factors, and antimicrobial resistance genes. We also demonstrate the potential of nanopore sequencing assemblies for rapid preliminary phylogenetic inference. Nanopore sequencing provides additional advantages as very low capital investment and footprint, and shorter (10 hours library preparation and sequencing) turnaround time compared to other NGS technologies.

MinION sequencing to genotype US strains of infectious laryngotracheitis virus.

Over the last decade the US broiler industry has fought long-lasting outbreaks of infectious laryngotracheitis (ILTV). Previously, nine genotypes (I-IX) of ILTVs have been recognized using the polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLP) method with three viral alleles (gB, gM and UL47/gG). In this study, the genotyping system was simplified to six genotypes by amplicon sequencing and examining discriminating single nucleotide polymorphisms (SNPs) within these open reading frames. Using phylogenomic analysis of 27 full genomes of ILTV, a single allele (ORF A/ORF B) was identified containing SNPs that could differentiate ILTVs into genotypes congruent with the phylogenetic partitioning. The allelic variations allowed for the cataloging of the 27 strains into 5 genotypes: vaccinal TCO, vaccinal CEO, virulent CEO-like, virulent US and virulent US backyard flocks from 1980 to 1990, correlating with the PCR-RFLP genotypes I/ II/ III (TCO), IV (CEO), V (virulent CEO-like), VI (virulent US) and VII/VIII/IX (virulent US backyard flock isolates). With the unique capabilities of third generation sequencing, we investigated the application of Oxford Nanopore MinION technology for rapid sequencing of the amplicons generated in the single-allele assay. This technology was an improvement over Sanger-based sequencing of the single allele amplicons due to a booster amplification step in the MinION sequencing protocol. Overall, there was a 90% correlation between the genotyping results of the single-allele assay and the multi-allele assay. Surveillance of emerging ILTV strains could greatly benefit from real-time amplicon sequencing using the single-allele assay and MinION sequencing. RESEARCH HIGHLIGHTS A multi-allelic assay identified nine ILTV genotypes circulating in the US Single-allele genotyping is congruent with whole genome phylogenetic partitioning US ILTV strains can be grouped into five genotypes using the single-allele assay The single-allele assay can be done using MinION sequencing of barcoded amplicons.

Proteome-wide Analysis of Protein Thermal Stability in the Model Higher Plant Arabidopsis thaliana.

Modern tandem MS-based sequencing technologies allow for the parallel measurement of concentration and covalent modifications for proteins within a complex sample. Recently, this capability has been extended to probe a proteome's three-dimensional structure and conformational state by determining the thermal denaturation profile of thousands of proteins simultaneously. Although many animals and their resident microbes exist under a relatively narrow, regulated physiological temperature range, plants take on the often widely ranging temperature of their surroundings, possibly influencing the evolution of protein thermal stability. In this report we present the first in-depth look at the thermal proteome of a plant species, the model organism Arabidopsis thaliana By profiling the melting curves of over 1700 Arabidopsis proteins using six biological replicates, we have observed significant correlation between protein thermostability and several known protein characteristics, including molecular weight and the composition ratio of charged to polar amino acids. We also report on a divergence of the thermostability of the core and regulatory domains of the plant 26S proteasome that may reflect a unique property of the way protein turnover is regulated during temperature stress. Lastly, the highly replicated database of Arabidopsis melting temperatures reported herein provides baseline data on the variability of protein behavior in the assay. Unfolding behavior and experiment-to-experiment variability were observed to be protein-specific traits, and thus this data can serve to inform the design and interpretation of future targeted assays to probe the conformational status of proteins from plants exposed to different chemical, environmental and genetic challenges.

Rapid virulence prediction and identification of Newcastle disease virus genotypes using third-generation sequencing.

BACKGROUND: Newcastle disease (ND) outbreaks are global challenges to the poultry industry. Effective management requires rapid identification and virulence prediction of the circulating Newcastle disease viruses (NDV), the causative agent of ND. However, these diagnostics are hindered by the genetic diversity and rapid evolution of NDVs. METHODS: An amplicon sequencing (AmpSeq) workflow for virulence and genotype prediction of NDV samples using a third-generation, real-time DNA sequencing platform is described here. 1D MinION sequencing of barcoded NDV amplicons was performed using 33 egg-grown isolates, (15 NDV genotypes), and 15 clinical swab samples collected from field outbreaks. Assembly-based data analysis was performed in a customized, Galaxy-based AmpSeq workflow. MinION-based results were compared to previously published sequences and to sequences obtained using a previously published Illumina MiSeq workflow. RESULTS: For all egg-grown isolates, NDV was detected and virulence and genotype were accurately predicted. For clinical samples, NDV was detected in ten of eleven NDV samples. Six of the clinical samples contained two mixed genotypes as determined by MiSeq, of which the MinION method detected both genotypes in four samples. Additionally, testing a dilution series of one NDV isolate resulted in NDV detection in a dilution as low as 101 50% egg infectious dose per milliliter. This was accomplished in as little as 7 min of sequencing time, with a 98.37% sequence identity compared to the expected consensus obtained by MiSeq. CONCLUSION: The depth of sequencing, fast sequencing capabilities, accuracy of the consensus sequences, and the low cost of multiplexing allowed for effective virulence prediction and genotype identification of NDVs currently circulating worldwide. The sensitivity of this protocol was preliminary tested using only one genotype. After more extensive evaluation of the sensitivity and specificity, this protocol will likely be applicable to the detection and characterization of NDV.

Draft Genome Sequences of Three Ochrobactrum spp. Isolated from Different Avian Hosts in Pakistan.

Here, we present the draft genome sequences of three Ochrobactrum sp. strains with multidrug-resistant properties, isolated in 2015 from a pigeon and two chickens in Pakistan.

Draft Genome Sequences of Five Novel Ochrobactrum spp. Isolated from Different Avian Hosts in Nigeria.

Here, we present the draft genome sequences of five multidrug-resistant novel Ochrobactrum species strains isolated from a pigeon, a duck, and chickens from Nigeria in 2009.

A robust and cost-effective approach to sequence and analyze complete genomes of small RNA viruses.

BACKGROUND: Next-generation sequencing (NGS) allows ultra-deep sequencing of nucleic acids. The use of sequence-independent amplification of viral nucleic acids without utilization of target-specific primers provides advantages over traditional sequencing methods and allows detection of unsuspected variants and co-infecting agents. However, NGS is not widely used for small RNA viruses because of incorrectly perceived cost estimates and inefficient utilization of freely available bioinformatics tools. METHODS: In this study, we have utilized NGS-based random sequencing of total RNA combined with barcode multiplexing of libraries to quickly, effectively and simultaneously characterize the genomic sequences of multiple avian paramyxoviruses. Thirty libraries were prepared from diagnostic samples amplified in allantoic fluids and their total RNAs were sequenced in a single flow cell on an Illumina MiSeq instrument. After digital normalization, data were assembled using the MIRA assembler within a customized workflow on the Galaxy platform. RESULTS: Twenty-eight avian paramyxovirus 1 (APMV-1), one APMV-13, four avian influenza and two infectious bronchitis virus complete or nearly complete genome sequences were obtained from the single run. The 29 avian paramyxovirus genomes displayed 99.6% mean coverage based on bases with Phred quality scores of 30 or more. The lower and upper quartiles of sample median depth per position for those 29 samples were 2984 and 6894, respectively, indicating coverage across samples sufficient for deep variant analysis. Sample processing and library preparation took approximately 25-30 h, the sequencing run took 39 h, and processing through the Galaxy workflow took approximately 2-3 h. The cost of all steps, excluding labor, was estimated to be 106 USD per sample. CONCLUSIONS: This work describes an efficient multiplexing NGS approach, a detailed analysis workflow, and customized tools for the characterization of the genomes of RNA viruses. The combination of multiplexing NGS technology with the Galaxy workflow platform resulted in a fast, user-friendly, and cost-efficient protocol for the simultaneous characterization of multiple full-length viral genomes. Twenty-nine full-length or near-full-length APMV genomes with a high median depth were successfully sequenced out of 30 samples. The applied de novo assembly approach also allowed identification of mixed viral populations in some of the samples.

A proteomic atlas of the legume Medicago truncatula and its nitrogen-fixing endosymbiont Sinorhizobium meliloti.

Legumes are essential components of agricultural systems because they enrich the soil in nitrogen and require little environmentally deleterious fertilizers. A complex symbiotic association between legumes and nitrogen-fixing soil bacteria called rhizobia culminates in the development of root nodules, where rhizobia fix atmospheric nitrogen and transfer it to their plant host. Here we describe a quantitative proteomic atlas of the model legume Medicago truncatula and its rhizobial symbiont Sinorhizobium meliloti, which includes more than 23,000 proteins, 20,000 phosphorylation sites, and 700 lysine acetylation sites. Our analysis provides insight into mechanisms regulating symbiosis. We identify a calmodulin-binding protein as a key regulator in the host and assign putative roles and targets to host factors (bioactive peptides) that control gene expression in the symbiont. Further mining of this proteomic resource may enable engineering of crops and their microbial partners to increase agricultural productivity and sustainability.